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ABSTRACT 

By searching through more than 10 satelhte-years of THEMIS and Cluster 
data, three rehable examples of parallel electric field turbulence in the undis- 
turbed solar wind have been found. The perpendicular and parallel electric field 
spectra in these examples have similar shapes and amplitudes, even at large scales 
(frequencies below the ion gyroscale) where Alfvenic turbulence with no paral- 
lel electric field component is thought to dominate. The spectra of the parallel 
electric field fluctuations are power laws with exponents near -5/3 below the ion 
scales (~0.1 Hz), and with a flattening of the spectrum in the vicinity of this 
frequency. At small scales (above a few Hz), the spectra are steeper than -5/3 
with values in the range of -2.1 to -2.8. These steeper slopes are consistent with 
expectations for kinetic Alfven turbulence, although their amplitude relative to 
the perpendicular fluctuations is larger than expected. 

Subject headings: plasmas — solar wind — turbulence 
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The magnetic field, velocity and density fluctuations in the so lar wind have been mea- 
sured for many decades ( lColemanlll968l : llntriligator &: Wolfelll970f ) and the majority of our 
knowledge about solar wind turbulence has come from measurer nents of these qu antities. 
Electric field fluctuations have only been used relatively recently ( jBale et al.l 120051 ) but can 
provide important additional information about its nature, such as the phase speed of the 
fluctuations and the their possible dispersion properties. In this Letter, we present the first 
measurement of the electric field fluctuations of solar wind turbulence in the direction parallel 
to the magnetic field, and discuss the implications of the measurements. 

At spacecraft-measured frequencies smaller than around 0.5 Hz, the power spectra of 
magnetic field, veloci ty and density fluctuations haye a power law form with an exponent 



around -1.5 to -1.7 f lMatthaeus fc GoldsteinI Il982l : Marsch fc Tul llQQOl : ISmith et al.l 12006 
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Chen et al.ll201ll : iBoldyrev et al.ll201ll ). This is usually interpreted as inertial range turbu- 
lence, in which there is a cascade of energy from large to small scales. At higher frequencies, 
corresponding to convected scales smaller than the ion gyroscale and inertial length, the 



magnetic field and den sity spectra are seen to steepen (IDenskat et al. 



19981 : IChen et al.l 120121 ) , which is often at tributed to dampi r ig ([Goldstein et al. 



further cascade of dispersive fluctuations (JGhosh et al.lll996l : ISchekochihin et al. 



19831: 



Leamon et al. 



199^ 



or a 



20091) 



The first measurement of the electric field spectrum in the spacecraft frame shows that 
in the inertial range, it has a -5/3 scaling, the same as the magnetic field, and at high 
frequencies, aroun d the ion gyroscale, may become enhanced compared to the magnetic field 
( iBale et al.l l2005l ). The ratio of electric to magnetic field fluctuations was used to infer 
the phase spee d, which was found to be consisten t with kinetic Alfven waves rather than 



whistler waves ( Bale et al. 20051: Saleni et a 



i n the magnetosheath (ISundkvist et al.ll2007l ). aurora (IChaston et a 



20121 ). S imilar measurements have been made 



( jEastwood et al.ll2009l ). Sim ulations of both gyrokinetic turbulence ( 



and Hall MHD turbulence (IDmitruk fc Matthaeud l2006l : iMatthaeus et al. 



20081) and magnetotai l 



Howes et al. 



2008 



2008, 



20111 ) 



20101 ) also 



show an enhancement in electric field fluctuations at ion kinetic scales, although comparison 
to the solar wind observations must be made carefully since, unlike the observations, the 
simulations are in a zero mean velocity frame. 

It has been noted that it is impor tant to consider th e measurement frame when in- 
terpreting electric field measurements (JKellogg et al.l l2006l ). since the electric field is not 
Galilean invariant. It has been shown that in the inertial range, the electric field in the 
spacecraft frame follows the same scaling as the magnetic field due to the fast convection 
past the spacecraft of the magnetic fluctuations ( IChen et al.ll201ll ). In the frame of zero 
mean velocity, which is perhaps more interesting from a theoretical viewpoint, the electric 
field takes a shallower scaling of -1.4. 

For scales smaller than the ion kinetic scales, the perpendicular electric field spectrum in 

—1/3 

the zero mean velocity frame is predicted to have a scaling oi kj_ , where k \ is the wavevec- 



tor perpendicular to the magnetic field direction (ISchekochihin et al.ll2009|). A more recent 
model, which includes strong intermittency, predicts a scaling of A;^ ( Boldyrev fc Perez 
20121 ). 



At sc ales larger than the ion kinetic scales, th e turbulence is thought to be predominantly 
Alfvenic ( iBelcher fc Davislll97ll : lBale et al.ll2005l ). with no electric field fluctuations parallel 
to the magnetic field, although such fluctuations may arise from the compressive component 
of the turbulence. At smaller scales, i.e. frequencies greater than 0.5 Hz the Alfvenic 



turbulence is th ought to develop parallel electric field fluctuations ( iBian fc Kontarl 12010 



Bian et al.ll2010l ). Understanding this parallel electric field is important for several reasons. 
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At the ion and electron kinetic scales it is thought that the turbulence is damped, which may 
occur partly by Landau darn ping, in which the particles interact with the parallel electric field 



fiuctuations (lQuataertlll998f). These fluctu ations may also be important for the non-thermal 



acceleration of particles (iBian et al.ll2010l . and references therein). Finally, measuring these 
fluctuations gives us a new way to probe the fundamental physics of kinetic scale plasma 
turbulence. 



2. Data 

Criteria for a successful parallel electric fleld spectral measurement in the quiet solar 
wind are stringent. A flrst requirement is that data be transmitted at a high data rate 
(8192 samples/second for THEMIS) and high gain in order to obtain spectral information 
at frequencies above a few Hz. This high data rate transmission occurred less than 0.1% 
of the time. A second criterion is that, during such high rate transmissions, the spacecraft 
cannot be in the tail, the foreshock, or the lunar wake. For THEMIS-C in lunar orbit, 
this happened about 20% of the time. A next requirement is that the magnetic fleld vector 
must be essentially in the spacecraft spin plane (which is approximately the ecliptic plane 
for THEMIS and Cluster). This is because electric fleld measurements along the spin axis 
were not made on Cluster and were not of sufficient quality on THEMIS due to its short on- 
axis antennas. The magnetic field was in the spin plane infrequently, perhaps a few percent 
of the time. A last requirement is that the spin plane magnetic field be perpendicular to 
the Sun-Earth line (in approximately the ecliptic Y-direction) because the sunward electric 
field component is perturbed by each of the four sensors rotating into the wake of the solar 
wind fiow every spin period. The combination of all these requirements led to finding two 
THEMIS examples and one Cluster example in searching about three years of data from the 
four Cluster and five THEMIS spacecraft. This data is presented below. 

Parallel electric field spectra were obtained in two different ways with the requirement 
of a good measurement being that the two spectral estimates agreed. These two estimates 
were made in two different coordinate systems, the despun spacecraft coordinates (dsc) and 
the field-aligned coordinates (fac) . The dsc coordinates are aligned with spacecraft axes such 
that Zdsc is along the spin axis (approximately the ecliptic normal) , Xdsc is in the spin plane 
and pointing sunward (approximately in the X-direction in GSE coordinates), and Ydsc, also 
in the spin plane, pointing generally westward. Because the Ydsc direction is the direction 
of the magnetic field in the three successful events, the parallel electric field in these events 
was also in the Ydsc direction. By contrast, in the field- aligned coordinate system, Zfac 
is the direction parallel to the instantaneous magnetic field, Xfac is perpendicular to the 
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instantaneous magnetic field, pointing generally sunward, and Yfac is the third direction 
in a right-handed coordinate system. Thus, on average, when the magnetic field is nearly 
in the Ydsc direction, EYdsc should be equal to EZfac, but they may be unequal in detail 
for the following reasons. The field-aligned coordinate system fiuctuates in space because 
the magnetic field direction varies. These angular fiuctuations were less than one or a few 
degrees for the examples shown. On the other hand, EYdsc is in a coordinate system 
fixed in space. Even so, EYdsc may contain components of the perpendicular electric field 
because the average and instantaneous directions of the magnetic field can be different. That 
these differences between EYdsc and EZfac are small and that either measurement gives the 
parallel electric field spectrum will be shown in the figures that follow. It is noted that a 
Galilean transformation to the plasma frame has not been made because the parallel electric 
field is Galilean invariant. 



Fig. la shows an example of magnetic field measurements (jAuster et al.l l2008l ) on the 
THEMIS-C spacecraft from 08:55:10 UT to 09:04:10 UT on April 21, 2012, when the space- 
craft was approximately on the Sun-Earth line and 60 Earth radii upstream and the plasma 
beta was 0.54. Because the magnetic field was in the Ydsc direction, this event satisfies 
the earlier selection criteria and the parallel electric field was in the Ydsc direction. Any 
leakage of the perpendicular electric field into the Ydsc component can be at most a few 
percent because the ratio of BXdsc or BZdsc to BYdsc in Fig. la is a few percent. (We 
shall see below that the parallel and perpendicular spectra are about equal so the parallel 
electric field spectrum cannot be due to leakage of t he perpendicular electric field). Fig. lb 



shows EYdsc time-domain data (JBonnell et al.ll2008l ) at a rate of 128 samples/second during 



the ten minute interval. It is noted that the turbulence level during the first part of the 
data interval was larger than near the end of the interval. This is not unusual for such time 
domain data. The spectral slopes through the first and last halves of the interval were the 
same, so the data presented below is the average over the entire interval. 

Parallel (EYdsc and EZfac) and perpendicular (EXdsc) electric field spectra are shown 
in Fig. 2a for the time interval of Fig. 1. An important demonstration of the validity of 
these spectral measurements is that the two parallel spectra agree. It is also noted that the 
perpendicular power spectrum, which is plotted a factor of 100 lower on the vertical scale 
to separate it from the roughly equal power parallel spectrum, had a magnitude and shape 
similar to that of the parallel spectrum. Wake effects or other perturbations of the waveform 
of Fig. lb were less than about 3 millivolts/meter (mV/m) and they were removed by notch 
filtering the waveform of Fig. lb at the fundamental and 11 harmonics to remove 1.2% of 
the frequency space between 10~^ and 64 Hz. Even so, it is possible that the spectra outside 
the region of these narrow spikes were perturbed by the presence of the wake. To check 
this possibility, half of the total data in the time domain waveform was removed by deleting 
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regions around each spike, and the spectral shape was the same as when all the data was 
used. 

Within the 10 minute time interval whose spectra cover below 0.01 to 5 Hz, data was 
also fed to the ground through a higher gain, higher data rate, AC coupled channel that 
produced the spectra at frequencies of 5 to 1000 Hz in Fig. 2a. The spectra of the total 
magnetic field and the trac e of its correlatio n spectral matrix are given in Fig. 2b, with 



fluxgate (JAuster et al.l 120081 ) and search coil (JRoux et al.l l2008l ) instruments providing, re 



spectively, the lower and higher frequency portions of the data. The peak in the magnetic 
field spectra near 50 Hz is due to a physical process whose discussion is beyond the scope of 
this paper. This peak is also seen at times in the higher frequency electric field spectra and 
a time interval was chosen for analysis that did not contain this spectral peak. Because the 
higher frequency data collection was triggered by an on-board algorithm that selected the 
largest amplitude events, the higher frequency spectra contained more power than their lower 
frequency counterparts. To compensate for this effect, the powers in these higher frequency 
components were decreased in order that the high and low frequency data fit together. It is 
important to note that the power in the perpendicular electric field is about equal to that in 
the parallel electric field, which shows that the observed parallel electric field does not result 
from contamination of the perpendicular electric field. 

The power law exponents of different segments of the electric and magnetic field spectra 
are given near each segment of each plot in Fig. 2. With the exception of the lowest 
frequency portion of the perpendicular electric field spectrum (which may differ from the 
other spectra due to inadequate data coverage), the low frequency spectra all have power law 
slopes of about -5/3. This includes the parallel electric field spectra that, for pure Alfvenic 
turbulence, would have no power at these lowest frequencies. All the spectra are relatively 
flat near the proton gyrofrequency (the vertical dotted curve) and have spectral slopes of 
about -5/3 above this frequency until, above a few Hz, the spectra become steeper and have 
power law slopes between about -1.9 and -2.8. The three vertical dotted lines are, from 
lowest to highest frequency, the proton gyrofrequency, the proton gyroscale (the solar wind 
velocity divided by 27r times the proton gyroradius), and the electron gyroscale. 

A second example of electric and magnetic field spectra is given in Fig. 3 for THEMIS- 
C data collected on October 14, 2012 from 10:27 UT to 10:38 UT when the plasma beta 
was 1.1. The features of these plots are the same as those described for the data of Fig. 2. 
Namely, the two estimates of the parallel electric field power agree, the parallel spectrum has 
somewhat more power than the perpendicular spectrum, and the three power spectra have 
the same shape. It is emphasized again that the parallel electric field power at the lowest 
frequencies was very much non-zero. Also, the parallel electric field, perpendicular electric 
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field, magnetic field trace, and total magnetic field all have power law slopes with spectral 
index about -5/3 at the lowest frequencies. The spectra are all flatter near the ion gyroscale 
and have slopes of about -5/3 above this frequency until, above a frequency of a few Hz, the 
slopes become steeper to lie between about -2.1 and -2.5. 



Spectral data from Cluster-3 electric (JGustafsson et al.lll997l ) and magnetic field ( JBalogh et al 



19971 ) measurements covering 11.5 minutes at a distance of 19 Earth radii, a magnetic lat- 
itude of 19 degrees, and a magnetic local time of 10:15, are given in Fig. 3. For this case 
also, the magnetic field was in the Ydsc direction and the two parallel electric field spectra 
in Fig. 3a are in essential agreement. Because the antennas on Cluster were 88 meters, 
tip-to-tip, wake perturbations were much smaller than those observed on THEMIS whose 
pair of spin plane antennas were 40 and 50 meters, tip-to-tip. Nevertheless, precautions 
taken in treating wake effects in the THEMIS data were also used in analyzing the Cluster 
event. Cluster was in the undisturbed solar wind during this time interval, as expected from 
the magnetic field being in the Y-direction and as shown by the lack of foreshock plasma 
in the electrostatic analyzers and the presence of Langmuir turbulence in the Whisper in- 
strument. The frequency coverage for this event extended to only 1 Hz because there was 
no higher frequency data collection. For this event, the parallel electric field power was an 
order-of-magnitude larger than the perpendicular electric field power. Otherwise the parallel 
and perpendicular spectra were essentially identical and they had power law slopes over the 
full frequency range of about -5/3 except in the vicinity of the ion gyrofrequency where the 
slopes were flatter. 



3. Discussion 

For the three examples of parallel electric field measurements in the solar wind, the 
parallel electric field estimated as EYdsc agreed with that estimated as EZfac, thereby 
vahdating the methods used to obtain the parallel electric field turbulence spectra. In general, 
the parallel electric field powers were comparable to or greater than the perpendicular powers, 
which rules out the possibility that the parallel fields were due to contamination from the 
perpendicular fields. Both the electric and magnetic field spectra exhibited plateaus near the 
ion gyrofrequency (~0.1 Hz) and had exponential slopes of about -5/3 below and above this 
frequency. The two events having spectra that extended to higher frequencies had spectra 
above a few Hz that became steeper with slopes between -1.9 and -2.8. 

For low frequencies, below 0.1 Hz, the slope of the parallel electric field sp ectrum matches 



that of the magnetic fie ld, which is not expected for pure Alfvenic turbulence (JBian fc Kontar 



2010l : iBian et al.ll2010l ). This may be due to the parallel electric field fluctuations from the 
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slow mode like com ponent of the turbulence, which contains around 10% of the energy 
( JHowes et al.l l2012l and references therein) or could be due to the non-linear turbulence 
differing from the predi ctions of the linear wave modes. Similarly to the density spectrum 
( IChen et al.ll2012l |2013[ ). the parallel electric field spectrum also displays a local flattening 
around the ion sc ales. It is possible that this is for a similar reason to that discussed by 
Chen et al.l ( l2013l ): the parallel electric field fluctuations from the kinetic Alfven turbulence 
take over from the slow mode at this scale, causing the observed plateau. Smaller plateaus 
are also seen in the trace magnetic field spectra; the reason for this is not known, since 
they are not usually seen, although this feature may be related to the short interval length, 
meaning that the low frequencies are not well measured. 

The parallel elect r ic field spectrum at high frequencies is steeper than some predictions 
( iBian &: Kontarl l2010t iBian et al.l l2010f ) although it is in a similar range to some recent 
simulations of strong kinetic Alfven turbulence (S. Boldyrev, private communication). Its 
amplitude relative to the perpendicular spectrum, however, is larger than expected. This 
may be due to other types of fluctuations being present at these scales, or to the non-linear 
turbulence differing from the linear wave mode predictions. The presence of the significant 
parallel electric field fluctuations has implications for how the turbulence heats the solar 
wind; in particular, it means that in addition to transit-time dampi ng from the pa rallel 
magnetic field fluctuations, there may be significant Landau damping (lQuataertlll998[ ). 



In the future it would be desirable to obtain a larger number of intervals to check 
whether the results presented here are valid at a statistical level. Due to the lack of current 
data that meets the stringent conditions required for the analysis, a dedicated solar wind 
turbulence mission may be required for this. 

This work was supported by NASA contracts NAS5-02099-09/12, NNN06AA01C and 
NASA grants NNX13AE24G, NNX09AE41G-1/14. 



REFERENCES 

Auster, H. U., Glassmeier, K. H., Magnes, W., et al. 2008, Space Sci. Rev., 141, 235 

Bale, S. D., Kellogg, R J., Mozer, F. S., Horbury, T. S., & Reme, H. 2005, Phys. Rev. Lett., 
94, 215002 

Balogh, A., Dunlop, M. W., Cowley, S. W. H., et al. 1997, Space Sci. Rev., 79, 65 

Belcher, J. W., & Davis, L. 1971, J. Geophys. Res., 76, 3534 



Bian, N. H., & Kontar, E. P. 2010, Phys. Plas., 17, 062308 

Bian, N. H., Kontar, E. P., & Brown, J. C. 2010, A&A, 519, AIM 

Boldyrev, S., & Perez, J. C. 2012, ApJ, 758, L44 

Boldyrev, S., Perez, J. C, Borovsky, J. E., & Podesta, J. J. 2011, ApJ, 741, L19 

Bonnell, J. W., Mozer, F. S., Delory, G. T., et al. 2008, Space Sci. Rev., 141, 303 

Chaston, C. C, Salem, C, Bonnell, J. W., et al. 2008, Phys. Rev. Lett., 100, 175003 

Chen, C. H. K., Bale, S. D., Salem, C, & Mozer, F. S. 2011, ApJ, 737, L41 

Chen, C. H. K., Howes, G. G., Bonnell, J. W., et al. 2013, Solar Wind 13 Proc. (in press), 
li?Xrv:1210.0127l 

Chen, C. H. K., Salem, C. S., Bonnell, J. W., Mozer, F. S., & Bale, S. D. 2012, 
Phys. Rev. Lett., 109, 035001 

Coleman, P. J. 1968, ApJ, 153, 371 

Denskat, K. U., Beinroth, H. J., & Neubauer, F. M. 1983, J. Geophys., 54, 60 

Dmitruk, P., & Matthaeus, W. H. 2006, Phys. Plas., 13, 042307 

Eastwood, J. P., Phan, T. D., Bale, S. D., & Tjulin, A. 2009, Phys. Rev. Lett., 102, 035001 

Ghosh, S., Siregar, E., Roberts, D. A., & Goldstein, M. L. 1996, J. Geophys. Res., 101, 2493 

Goldstein, M. L., Roberts, D. A., & Fitch, C. A. 1994, J. Geophys. Res., 99, 11519 

Gustafsson, G., Bostrom, R., Holback, B., et al. 1997, Space Sci. Rev., 79, 137 

Howes, G. G., Bale, S. D., Klein, K. G., et al. 2012, ApJ, 753, L19 

Howes, G. G., Borland, W., Cowley, S. C, et al. 2008, Phys. Rev. Lett., 100, 065004 

Howes, G. G., TenBarge, J. M., Borland, W., et al. 2011, Phys. Rev. Lett., 107, 035004 

Intriligator, B. S., & Wolfe, J. H. 1970, ApJ, 162, L187 

Kellogg, P. J., Bale, S. B., Mozer, F. S., Horbury, T. S., & Reme, H. 2006, ApJ, 645, 704 

Leamon, R. J., Smith, C. W., Ness, N. F., Matthaeus, W. H., & Wong, H. K. 1998, J. Geo- 
phys. Res., 103, 4775 



-9- 

Marsch, E., & Tu, C.-Y. 1990, J. Geophys. Res., 95, 11945 

Matthaeus, W. H., & Goldstein, M. L. 1982, J. Geophys. Res., 87, 6011 

Matthaeus, W. H., Servidio, S., & Dmitruk, R 2008, Phys. Rev. Lett., 101, 149501 

Matthaeus, W. H., Servidio, S., & Dmitruk, R 2010, AIP Conf. Proc, 1216, 184 

Quataert, E. 1998, ApJ, 500, 978 

Roux, A., Le Contel, O., Coillot, C, et al. 2008, Space Sci. Rev., 141, 265 

Salem, C. S., Howes, G. G., Sundkvist, D., et al. 2012, ApJ, 745, L9 

Schekochihin, A. A., Cowley, S. C., Borland, W., et al. 2009, ApJS, 182, 310 

Smith, C. W., Hamilton, K., Vasquez, B. J., & Leamon, R. J. 2006, ApJ, 645, L85 

Sundkvist, B., Retino, A., Vaivads, A., & Bale, S. B. 2007, Phys. Rev. Lett., 99, 025004 



This preprint was prepared with the AAS E^TJtjX macros v5.2. 



-10- 




Eydsc 
mV/m 



-1 







1 1 



3 6 

TIME AFTER 0855:10 UT, minutes 



Fig. 1. — Three components of the magnetic field (Fig. la) and the Y-component of the 
electric field (Fig. lb) during THEMIS-C spacecraft observations from 08:55:10 UT to 
09:04:10 UT on April 21, 2012, when the spacecraft was approximately on the Sun-Earth 
line and 60 Earth radii upstream. Because the magnetic field was in the Y-direction, the 
Y-component of the electric field was the parallel component of the electric field. 
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Fig. 2. — Spectra of the parallel electric field measured two different ways, the perpendicular 
electric field, the magnetic field trace and the total magnetic field during the data period 
illustrated in Fig. 1. The validity of the measurement is shown by the fact that the two 
parallel electric field spectra agree. In addition, the parallel and perpendicular electric field 
spectra have the same magnitudes, showing that the parallel electric field did not arise from 
contamination from the perpendicular electric field. The power law slopes (the numbers 
near each segment of the plots) are in agreement with expectations except that the parallel 
electric field spectral amplitude is unexpectedly large. The three vertical dotted hues are, 
from lowest to highest frequency, the proton gyrofrequency, the proton gyroscale, and the 
electron gyroscale. 
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Fig. 3. — Spectra similar to those in Figure 2 for data collected on THEMIS-C in the 
undisturbed solar wind on October 14, 2012 from 10:27 to 10:38 UT. The conclusions from 
Figure 2 are the same as those found from this figure. 
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Fig. 4. — Spectra similar to those in Figures 2 and 3 for data collected in the undisturbed 
solar wind by the Cluster-3 satellite on February 10, 2001 06:17:30 to 06:29 UT. The con- 
clusions from Figures 2 and 3 are the same as those found from this figure. 



